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Abstract 

We  have  created  a  timescale  that  joins  the  short-term  stability  of  several  hydrogen  masers 
with  the  long-term  capabilities  of  an  ensemble  of  cesium  frequency  standards.  Control  theory  is 
utilized  in  a  system  design  that  combines  frequency  standards  with  varied  properties.  The 
system  steers  a  maser  ensemble  with  respect  to  a  cesium  ensemble  while  minimally  perturbing 
the  maser  short-term  performance.  Results  are  given  from  both  simulated  and  measured  data. 
Systems  are  designed  using  linear  quadratic  Gaussian  (LQG)  and  pole-placement  techniques. 


INTRODUCTION 

One  of  the  many  challenges  in  the  construction  of  a  timescale  (mean)  is  the  combination  of  frequency 
standards  with  dissimilar  stochastic  properties.  A  common  example  for  timing  institutions,  and  one  that 
is  addressed  in  this  paper,  is  that  of  combining  hydrogen  masers  and  cesium  standards.  The  hydrogen 
masers  are  approximately  an  order  of  magnitude  more  stable  than  the  cesium  standards  in  the  short  term, 
while  the  cesium  standards  show  better  long-term  characteristics.  There  are  roughly  60  cesium  standards 
and  15  hydrogen  masers  contributing  to  the  U.S.  Naval  Observatory  (USNO)  operational  timescale  at  any 
given  time.  These  standards  are  spread  throughout  several  locations  at  USNO  in  Washington,  DC. 


BASIC  DESIGN 

Figure  1  is  a  block  diagram  that  shows  the  overall  design  of  the  system.  It  is  comprised  of  two  coupled 
subsystems.  One  subsystem  is  created  by  steering  a  maser  mean  (MM)  to  a  steered  cesium  mean  (CM) 
[11-  The  cesium  mean  steers  are  those  used  to  align  it  with  UTC. 

The  other  subsystem  creates  the  physical  output  by  steering  the  output  of  a  maser  to  the  MM  with  a 
precision  frequency  synthesizer  referred  to  as  an  AOG. 

The  filter  utilized  in  both  subsystems  is  a  Kalman  filter  [2-4]. 

The  present  Master  Clock  (MC)  control  design  steers  a  maser-based  AOG  signal  once  a  day  to  a 
dynamically  weighted  timescale  [1].  The  prototype  designs  use  hourly  steers  and  the  masers  have  had 
deterministic  rates  and  drifts  with  respect  to  the  cesium  mean  removed. 
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Figure  1.  Design  Block  Diagram. 


STATE-SPACE  MODEL 

The  performance  of  frequency  standards  may  be  described  using  state-space  modeling  techniques  [3,4,5]. 
The  state-space  model  for  the  system  described  in  Figure  1  is 

Xt+T=<DXt+But 
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x  is  phase  difference,  v  is  frequency  difference,  and  the  control  ut  =  -GXt ,  where  G  is  the  control  gain. 
The  above  equation  can  now  be  rewritten  as  Xt+T  =  (<D  -  BG)  Xt . 


PROTOTYPE  SYSTEM  DESIGN 

The  linear  quadratic  Gaussian  (LQG)  design  technique  was  utilized  for  the  first  implementation  of  this 
system.  In  LQG  theory  [1,3, 5-7],  the  system  is  modeled  with  a  linear  state-space  equation  with  Gaussian 
noise  characteristics.  The  optimal  control  minimizes  a  quadratic  cost  function: 

J  =  Z[X^WQXt  +  ^WruJ 

t 

where, 

x 

Xt  = 

[y. 


phase  offset 
frequency  offset 


and 


u  is  a  vector  representing  the  control  ut  =  -GXt . 
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The  cost  matrix  weighting  for  state  offsets  is  Wq  and  the  cost  for  the  control  effort  WR.  The  solution  to 
minimizing  the  quadratic  cost  function  J  can  be  shown  to  yield  [1,3, 5-7]: 

G  -  (BTKB  +  Wr)_1BtK®  , 
where  K  is  a  solution  to  the  steady-state  Riccati  equation 

K  =  +  WQ  -$tKB(BtKB  +  Wr)-1BtKG> 


Weighting  matrices  were  tested  using  simulations  until  a  control  system  was  obtained  that  provided 
minimal  short-term  steering  perturbations  along  with  small  time  and  frequency  offsets.  The  weighting 
matrices  selected  based  on  the  simulations  were 
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The  control  gain  calculated  from  the  Riccati  equation  is 
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G  = 
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Figures  2  and  3  show  the  performance  of  the  initial  design  as  measured  during  actual  real-time 
operational  system  testing.  The  time  and  frequency  offsets  were  acceptable,  but  the  frequency  stability  of 
the  system  for  averaging  times  greater  than  30,000  seconds  was  not  making  an  improvement  over  the 
present  operational  master  clock  steering  design. 


Figure  2.  Measured  Time  Difference  Data.  Figure  3.  Stability  Comparison. 
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SYSTEM  REDESIGN 


Although  the  system  performed  well,  it  was  decided  to  redesign  the  control  in  an  attempt  to  improve  the 
longer-term  frequency  stability  of  the  physical  output.  The  redesign  was  done  using  pole-placement 
techniques  in  order  to  attempt  to  get  a  better  handle  on  the  system  time  constants.  The  two-state  model 
for  a  frequency  standard  steered  utilizing  discrete  frequency  steps  is  given  by 


Xt+x  -  t&Xf  +  But  , 


where  X 
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The  eigenvalues  of  A  correspond  to  the  system  poles  and  are  related  to  the  response  of  the  system.  The 
relationship  between  the  control  gain  and  the  system  poles  can  be  solved  from 
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where  z  relates  to  the  z-transform  and  I  is  the  identity  matrix  [7].  By  equating  terms  in  the  above 
equation  with  the  eigenvalue  equation,  one  can  solve  for  the  control  gains. 

If  we  assume  each  subsystem  is  critically  damped,  then 

g*  =  r-10  -  exp  [-  r/Tc])2 
gy  =  1  “  exp  [-2  r/TJ, 

where  ris  the  constant  interval  between  steers  and  Tc  is  the  approximate  time  constant. 

We  based  the  control  gain  time  constants  for  the  coupled  system  on  the  assumption  that  it  was  actually 
separated  into  two  independent  two-state  systems.  The  time  constant  chosen  for  the  maser  mean  to 
cesium  mean  subsystem  was  30  days  and  the  time  constant  for  steering  the  AOG  to  the  maser  mean  was  4 
days.  The  control  gain  was  calculated  to  be 


G  - 


-  3.0e  -  8 
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-  0.02  0  0 
0  -  5.35e  - 10  -0.0027 


Figures  4-6  show  the  response  of  the  two  system  designs  to  initial  5  x  10'15  frequency  offsets  of  the 
cesium  mean,  maser  mean,  and  the  AOG. 
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Figure  4.  Responses  to  MM  frequency  offset. 
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Figure  5.  Responses  to  CM  frequency  offset. 
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Figure  6.  Responses  to  AOG  frequency  offset. 


A  comparison  of  the  measured  performance  of  the  initial  design  and  the  simulated  performance  of  the 
system  redesign  is  shown  in  Figure  7.  The  redesign  shows  improved  mid  to  long-term  stability 
performance.  Figure  8  shows  the  stability  of  the  underlying  maser  and  illustrates  a  decrease  in  the  control 
effort  of  the  simulated  redesign  compared  to  the  actual  steers  implemented  with  AOG  12. 


CONCLUSION 

Two  control  theory  design  techniques  were  utilized  in  a  scheme  for  combining  frequency  standards  that 
exhibit  different  frequency  stability  properties  to  make  a  physically  realizable  output  of  a  timescale.  The 
techniques  are  both  appropriate  for  designing  this  system.  The  first  prototype  design  gave  a  35% 
improvement  as  measured  by  the  Allan  deviation  at  2  days  compared  to  the  present  Master  Clock.  The 
redesign  promises  to  also  improve  the  mid-term  performance  without  the  long-term  degradation  seen  by 
the  prototype  design.  We  plan  to  continue  investigating  the  problem  and  hope  to  gain  more  insight  into 
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improving  our  understanding  of  how  system  response  characteristics  influence  the  overall  performance. 
This  will  likely  include  tests  on  how  controlling  each  maser  individually  affects  the  system. 


Figure  7.  Stability  comparison  of  different  designs.  Figure  8.  Control  effort  comparison. 
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